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Hadronic decays of the r lepton provide a clean environment to study hadron dynamics in an energy regime 
dominated by resonances. The interesting information is captured in the spectral functions. Recent results from 
ALEPH on exclusive channels are presented, with emphasis on the tttt'' final state which plays a crucial role for 
the determination of the hadronic contribution to the muon anomalous magnetic moment. A comparison between 
27r spectral functions obtained in r decays (after corrections for isospin- breaking) and e+e~ annihilation reveals 
some discrepancy in the line shape of the p resonance which can be attributed to different pole mass values for 
the charged and neutral p's, which are determined through a robust fitting procedure. However, after applying 
this correction, the normalization of the two spectral functions diff'er by 3.3%. Inclusive spectral functions are the 
basis for QCD analyses, which deliver an accurate determination of the strong coupling constant, quantitative 
information on nonperturbative contributions and a measurement of the mass of the strange quark. 



1. Introduction 

Hadrons produced in r decays are born out of 
the charged weak current, i.e. out of the QCD 
vacuum, which guarantees that hadronic physics 
factorizes. These processes are then completely 
characterized for each decay channel by a spec- 
tral function which can be directly extracted from 
the invariant mass spectra of the final state. Fur- 
thermore, the produced hadronic systems have 
7=1 and spin-parity = O"*",!" (V) and 
= 0^, 1+ (A). Isospin symmetry (CVC) con- 
nects the T and e~*'e~ annihilation spectral func- 
tions. 

Hadronic r decays are a clean probe of hadron 
dynamics in an interesting energy region dom- 
inated by resonances. However, perturbative 
QCD can be seriously considered due to the rel- 
atively large r mass. Samples of ~ 4 x 10^ mea- 
sured decays are available in each LEP experi- 
ment and CLEG. Conditions for low systematic 
uncertainties are particularly well met at LEP: 
measured samples have small non-r backgrounds 
(~ 1%) and large selection efficiency (92%), for 
example in ALEPH. 



2. Spectral Functions from ALEPH 

Preliminary spectral functions based on the 
full LEPl statistics are available from ALEPH. 
The corresponding results [2] for the branching 
fractions which provide the absolute normaliza- 
tion for the spectral functions are given in Ta- 
ble 1. The analysis uses an improved treatment 
of photons as compared to the published analyses 
based on a reduced sample [3,4]. Spectral func- 
tions are unfolded from the measured mass spec- 
tra after background subtraction using a mass- 
migration matrix obtained from the simulation 
in order to account for detector and reconstruc- 
tion biases. Backgrounds from non-rr events are 
small (< 1% and subtractions are dominated by 
r decay feedthroughs. 

The spectral functions are separated into vec- 
tor and axial- vector components according to the 
number of pions in the final state. As for final 
states involving KK pairs, specific input is re- 
quired to achieve the V — A separation [5] . 
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Branching fractions in r decays from the ALEPH 
experiment [2,5]. Apart from the two leptonic 
channels, the other modes are labeled according 
to their nonstrange vector (V), nonstrange axial- 
vector (A), and strange (S) hadronic final states. 
The u} decay modes marked (*) are electromag- 
netic (7r°7, TT+TT"). The branching fractions for 
7r~7r~7r"'"37r*^ and aj~(— *■ 7r~7) are estimates, while 
those for (3-77)^77 and (37r)^a; are from CLEO [6]. 
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Figure 1. The comparison between ALEPH and 
CLEO spectral functions in the tttt" mode. 



3. The 27r Vector State 

3.1. The Data 

The spectral function from r VtT^^tt''^ from 
the full-LEPl ALEPH analysis (- 10^ events) is 
available. It is in good agreement with the re- 
sults from CLEO [7] as shown in Fig. 1. The 
statistics is comparable in the two cases, however 
due to a flat acceptance in ALEPH and a strongly 
increasing one in CLEO, ALEPH data are more 
precise below the p peak, while CLEO is more 
precise above. Note that, due to the unfolding 
procedure, neighbouring data points are strongly 
correlated. 

3.2. TTTT Spectral Functions and tt Form Fac- 
tors 

It is useful to carefully write down all the fac- 
tors involved in the comparison of e^e~ and r 
spectral functions in order to make explicit the 
possible sources of CVC breaking. On the e+e" 
side we have 

cr(e"'"e~ — > 7r"'"7r~) = Vq{s) (1) 
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where /^g (s) is the threshold kinematic factor and 
F^(s) the pion form factor. On the r side, the 
physics is contained in the hadronic mass distri- 
bution through 
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where Vud = 0.9748 ± 0.0010 denotes the CKM 
weak mixing matrix element [1] and -Sew ac- 
counts for clectroweak radiative corrections (see 
below). SU(2) symmetry implies V-{s) = vq(s). 

Experiments on r decays measure the rate 
inclusive of radiative photons, i.e. for t — s- 
j>'-r7r7r°(7). The measured spectral function is thus 
v*_{s) = V-{s) G{s), where G{s) is a radiative 
correction, computed using scalar QED. 

Several levels of SU(2) breaking can be identi- 
fied: 

• electroweak radiative corrections to r decays 

are contained in the Sew factor [8,9] which 
is dominated by short-distance effects. As 
such it is expected to be weakly dependent 
on the specific hadronic final state, as veri- 
fied in the t — > (tt, K)i^t channels [10]. Re- 
cently, detailed calculations have been per- 
formed for the TTTT^ channel [11] which also 
confirm the relative smallness of the long- 
distance contributions. The total correction 
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where S^^ is the leading-log short-distance 
electroweak factor (which vanishes for lep- 
tons) and S^^'^^^ are the nonleading elec- 
tromagnetic corrections. The latter correc- 
tions are calculated in Ref. [9] at the quark 



level and in Ref. [11] at the hadron level for 
the tttt" decay mode, and in Ref. [8,9] for 
leptons. The total correction amomits [14] 
to S'f'{^" = 1.0198 ±0.0006 for the inclusive 
hadron decay rate and S^^r = (1.0232 ± 
0.0006) S'l;'^(s) for the tttt^ decay mode, 
where S]J}yj{s) is an s-dependent radiative 
correction [11]. 

• the pion mass splitting breaks isospin sym- 
metry in the spectral functions [12,13] since 
/3_(s)7^/3o(s). 

• symmetry is also broken in the pion form 
factor [12,11] from the tt mass splitting. 

• a similar effect is expected from the p mass 
splitting. The theoretical expectation [16] 
gives a limit (< 0.7 MeV), but it is re- 
ally only a rough estimate. Thus the ques- 
tion must be investigated experimentally, 
the best approach being the explicit com- 
parison of T and e+e~ 27r spectral func- 
tions, after correction for the other isospin- 
breaking effects. 

• explicit electromagnetic decays such as wy, 
777, aiid 77777 introduce some small dif- 
ference between the widths of the charged 
and neutral p's. 

• isospin violation in the strong amplitude 
through the mass difference between u and 
d quarks is expected to be negligible. 

3.3. Fitting the 27r Spectral Functions: the 

p mass splitting 

The 277 spectral function is dominated by 
the wide p resonance, parametrized following 
Gounaris-Sakurai [15] (GS) which takes into ac- 
count analyticity and unitarity properties. 

The pion form factor is fit with interfering am- 
plitudes from /9(770), /3'(1450) and p"(1700) vec- 
tor mesons with relative strengths 1, (3 and 7 (real 
numbers). A phase (pg is also considered, since 
the relative phase of the p and p' amplitude is a 
priori unknown. The much smaller relative ampli- 
tude for the p" is assumed to be real. In practice 
we fit F^{s) from e+e~ data and F~ (s) from the r 
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spectral function duly corrected for SU(2) break- 
ing, however only for the spectral function f3^ fac- 
tor, for the Sew factor, and for the long-distance 
radiative correction Gem{s). In this way, the 
mass and width of the dominating p resonance 
in the two isospin states can be unambiguously 
determined. All e+e" data are used, including 
the recently corrected precise results from CMD- 
2 [17,18]. On the t side, the accurate data from 
ALEPH and CLEO are used. 

The systematic uncertainties are included in 
the fits through appropriate covariance matrices. 
The p mass systematic uncertainty in r data is 
mostly from calibration (0.7 MeV for ALEPH 
and 0.9 MeV for CLEO). The corresponding un- 
certainties on the p width are 0.8 and 0.7 MeV. 
Both mass and width determinations are limited 
by systematic uncertainties in r data, however 
they can be safely assumed to be uncorrelated 
between ALEPH and CLEO. 

Due to the large event statistics, the fits are 
quite sensitive to the precise line shape and to 
the interference between the different amplitudes. 
Of course, p — uj interference is included for e+e^ 
data only and the corresponding amplitude {ap^,) 
is fit, together with its relative phase. The upper 
range of the fit is taken at 2.4 GeV^ in the t 
data in order to avoid the kinematic end-point of 
the T spectral function where large corrections are 
applied. The e'^e'~ data are fit up to 3.6 GeV^, so 
that the information on the p" comes essentially 
from e+e~ data. 

It is customary in the Novosibirsk analyses 
to define including the vacuum polarization 
(both leptonic and hadronic) in the photon prop- 
agator. Such a prescription is not desirable when 
studying the real hadronic structure of the pion, 
and incorrect when comparing to the r case where 
these contributions are absent. Removing vac- 
uum polarization in the fit to CMD-2 \F^\'^ data 
(Table 2) yields a p° mass 1.1 MeV lower. Sub- 
sequent fits are done excluding vacuum polariza- 
tion. 

It turns out that the resulting p masses and 
widths are quite sensitive to the strength of the 
p' and p" amplitudes, /3 and 7. So, depending 
on the type of fit, the derived values can ex- 
hibit some systematic shifts. The e+e~ and r 
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Table 2 

Fits to CMD-2 ]F^]2 data [17] using the Gounaris- 
Sakurai parametrization of the p, p' resonances 
for two definitions of the pion form factor: ex- 
cluding or including the vacuum polarization in 
the photon propagator. A mass shift of 1.1 MeV 
is observed for the p mass between the two cases. 
The values between squared brackets are fixed in 
the fits. The masses and widths are in MeV and 
(t)/3 ii^ degrees. 



fits yield significantly different p' amplitudes and 
phases: this is particularly true when a restricted 
energy range is used for the fit, as it is the case 
when only CMD-2 data are considered (only up 
to 960 MeV). In order to avoid this problem, a fit 
to both data sets is performed, keeping common 
values for the p' and p" parameters. In doing 
so, one neglects possible isospin-breaking effects 
for these states, which appears to be a reason- 
able assumption, as the p' and p" components 
can be considered as second-order with respect 
to the dominant investigated p amplitude. 

The differences between the masses and widths 
of the charged and neutral p's in the common fit 
are found to be 



(2.3 ± 0.8) MeV 
(0.1 ± 1.4) MeV 



(4) 
(5) 



The mass splitting is somewhat larger than the 
theoretical prediction (< 0.7 MeV) [16], but only 
at the 2 a level. The expected width splitting, 
not taking into account any p mass splitting, is 
(0.7 ± 0.3) MeV [11,14]. However, if the mass 
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difference is taken as an experimental fact, then a 
larger width difference would be expected. Prom 
the chiral model of the p resonance [19,11], one 
expects 

where ^Tem is the width difference from electro- 
magnetic decays (as discussed above), leading to a 
total width difference (2.1 ±0.5) MeV, marginally 
consistent with the observed value. 

Table 3 presents the results of the common fits, 
the quality of which can be inspected in Figs. 2 
and 3. 

Mass splitting for the p was in fact considered 
in our first paper where we proposed using pre- 
cise r data to compute hadronic vacuum polariza- 
tion integrals [12]. Using pre-CMD-2 e+e" data 
and ALEPH r, a combined fit was attempted 
which produced a mass splitting consitent with 
within an uncertainty of 1.1 McV. However, 
the form factor from e+e" data still contained 
the vacuum polarization contribution (1.1 MeV 
shift, as wc have seen) and we also discovered a 
normalization problem in our treatment of the 
r data in the combined fit. With the advent 
of precise CMD-2 data [17], it became appar- 
ent that differences were showing up between r 
and e+e~ form factors. A large part of the dis- 
crepancy was removed when CMD-2 re-analyzed 
their data [17]. Since the r results from ALEPH, 
CLEO and OPAL have been shown to be con- 
sistent within their quoted accuracy [14,18] and 
since preliminary results from the radiative return 
analysis of KLOE [20] are in excellent agreement 
with the corrected CMD-2 results, the question 
of the p mass splitting can be now more reliably 
investigated. Table 3 presents the results of the 
common fit, leaving m^- and mpo as free parame- 
ters, but fixing the relationship between F^- and 
FpO following Eq. (6). The results, which can be 
visualized in Figs. 2 and 3, do not provide a good 
description of the data with a probability of only 
0.6% to get a worse fit. 

Could a p mass splitting account for the 
discrepancy [14,18] between e+e~ and isospin- 
corrected r spectral functions? Unfortunately 
not! Correcting the r data for the mass shift us- 
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Table 3 

Combined fit to the pion form factor squared to 
T and e+e~ data, vacuum polarization excluded 
for the latter. The parametrization of the p, p' , 
p" line shapes follows the Gounaris-Sakurai form. 
For the p, only nip-, rUpO, and Fp- are fitted, 
while FpO is computed from Eq. (6). All mass 
and width values are in MeV and the phases are 
in degrees. The parameters related to p' and p" 
amplitudes are fitted, assuming they are identical 
in both data sets. 



ing the Gounaris-Sakurai parametrization with p, 

p' and p" amplitudes (but the results only de- 
pend on a Breit-Wigner-like behaviour for the 
p) improves the consistency of the e+e~ and r 
line shapes, but at the expense of an increased 
overall normalization difference; the corrected t 
data lie on average 3.3% above e+e~ hence the 
bad fit. The tau estimate of a^f-^'^'-' increases by 
5.4 10~^° [21], in larger disagreement with the 
e+e~ estimate [18]. 

4. Inclusive Nonstrange Spectral Func- 
tions 

The T nonstrange spectral functions have been 
measured by ALEPH [3,4] and OPAL [22]. The 
procedure requires a careful separation of vector 
(V) and axial-vector (A) states involving the re- 
construction of multi-7r° decays and the proper 
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Figure 2. The ratio of the e+e pion form fac- 
tor squared to |i^°(s)p from the combined fit in 
Table 3. 



Figure 3. The ratio of the t pion form factor 
squared to |i^~(s)p from the combined fit in Ta- 
ble 3. 



treatment of final states with a KK pair. The 
inclusive ALEPH V and A spectral functions are 
given in Figs. 4 and 5 with a breakdown of the 
respective contributions. The V + A spectral 
function, shown in Fig. 6 has a clear pattern con- 
verging toward a value above the parton level as 
expected in QCD. In fact, it displays a textbook 
example of global duality, since the resonance- 
dominated low-mass region shows an oscillatory 
behaviour around the asymptotic QCD expecta- 
tion, assumed to be valid in a local sense only for 
large masses. This feature will be quantitatively 
discussed in the next section. 

5. QCD Analysis of Nonstrange r Decays 

5.1. Motivation 

The total hadronic r width, properly normal- 
ized to the known leptonic width. 



r(r — > hadrons Vr) 



r(7 



(7) 



should be well predicted by QCD as it is an inclu- 
sive observable. Compared to the similar quan- 
tity defined in e+e" annihilation, it is even twice 
inclusive: not only are all produced hadronic 
states at a given mass summed over, but an inte- 
gration is performed over all the possible masses 
from to nir- 

This favourable situation could be spoiled by 
the fact that the scale is rather small, so that 
questions about the validity of a perturbative ap- 
proach can be raised. At least two levels are to be 
considered: the convergence of the perturbative 
expansion and the control of the nonperturbative 
contributions. Happy circumstances make these 
latter components indeed very small [23 25]. 

The perturbative expansion (FOPT) is known 
to third order [27]. A resummation of all known 
higher order logarithmic integrals improves the 
convergence of the perturbative series (contour- 
improved method FOPTci) [28]. As some ambi- 
guity persists, the results are given as an average 
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Figure 4. Preliminary ALEPH inclusive vector 
spectral function with its different contributions. 
The dashed line is the expectation from the naive 
parton model. 



Figure 5. Preliminary ALEPH inclusive axial- 
vector spectral function with its different contri- 
butions. The dashed line is the expectation from 
the naive parton model. 



of the two methods with the difference taken as 
a systematic uncertainty. The small nonpertur- 
bative contributions are parametrized using the 
Operator Product Expansion (OPE) [26]. 

5.2. Results of the Fits 

The QCD analysis of the r hadronic width has 
not yet been completed with the final ALEPH 
spectral functions. Results given below corre- 
spond to the published analyses with a smaller 
data set. 

The ratio Rr is obtained from measurements of 
the leptonic branching ratios: 

Rr = 3.647 ±0.014 (8) 

using a value B{t^ Ve^t) — (17.794 ± 

0.045)% which includes the improvement in ac- 
curacy provided by the universality assumption 
of leptonic currents together with the measure- 
ments of B(t~ — > e~ v^Vt), B{t^ —t ^.^ D^Vt) 
and the t lifetime. The nonstrange part of i?^ is 
obtained by subtracting out the measured strange 
contribution (see last section). The results of the 



fits are given in Table 4 for the ALEPH analysis. 

One notices a remarkable agreement within sta- 
tistical errors between the Q;s(m^) values using 
vector and axial-vector data. The total nonper- 
turbative power contribution to Rry+A is com- 
patible with zero within an uncertainty of 0.4%, 
that is much smaller than the error arising from 
the perturbative term. This cancellation of the 
nonperturbative terms increases the confidence 
on the as(TO^) determination from the inclusive 
iy -\- A) observables. 

The final result from ALEPH is : 

as{ml) = 0.334 ±0.007cxp± 0.021th (9) 

where the first error accounts for the experimen- 
tal uncertainty and the second gives the uncer- 
tainty of the theoretical prediction of R^ and the 
spectral moments as well as the ambiguity of the 
theoretical approaches employed. 
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Figure 6. The preliminary inclusive V + A 
nonstrange spectral function from ALEPH. The 
dashed line is the expectation from the naive par- 
ton model. 



5.3. Test of the Running of aa{s) at Low 
Energies 

Using the spectral functions, one can simulate 
the physics of a hypothetical r lepton with a mass 
y/so smaller than and hence further investi- 
gate QCD phenomena at low energies. Assum- 
ing quark-hadron duality, the evolution of Rr{so) 
provides a direct test of the running of as (so), 
governed by the RGE /3-function. On the other 
hand, it is a test of the validity of the OPE ap- 
proach in r decays. 

In practice, the experimental value for as{so) 
has been determined at every sq from the com- 
parison of data and theory. Good agreement is 
observed with the four-loop RGE evolution using 
three quark flavours (Fig. 7). The experimental 
fact that the nonperturbative contributions can- 
cel over the whole range 1.2 GeV^ < So < 
leads to confidence that the ag determination 
from the inclusive {V + A) data is robust. 
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0.330 ±0.014 ±0.018 
0.339 ±0.013 ±0.018 


0.020 ± 0.004 
-0.027 ± 0.004 


V+A 


0.334 ± 0.007 ± 0.021 


-0.003 ± 0.004 



Table 4 

Fit results of ^^(m^) and the OPE nonperturba- 
tive contributions from vector, axial-vector and 
{y ± A) combined fits using the corresponding 
ratios Rr and the spectral moments as input pa- 
rameters. The second error is given for theoretical 
uncertainty. 



5.4. Evolving as from m^) to A/z 

The evolution of the asirri^) measurement from 
the inclusive {V + A) observables based on the 
Runge-Kutta integration of the differential equa- 
tion of the renormalization group to N'^LO [29,31] 
yields for the ALEPH analysis 

as{Ml) = 0.1202±0.0008oxp±0.0024th±0.0010ovoi(10) 

where the last error stands for possible ambigu- 
ities in the evolution due to uncertainties in the 
matching scales of the quark thresholds [31]. The 
result (10) can be compared to the precise deter- 
mination from the measurement of the Z width, 
as obtained in the global electroweak fit. The 
variable Rz has similar advantages to Rt, but it 
differs concerning the convergence of the pertur- 
bative expansion because of the much larger scale. 
It turns out that this determination is dominated 
by experimental errors with very small theoretical 
uncertainties, i.e. the reverse of the situation en- 
countered in T decays. The most recent value [32] 
yields as{Mi) = 0.1183 ± 0.0027, in excellent 
agreement with (10). Fig. 7 illustrates well the 
agreement between the evolution of as(m^) pre- 
dicted by QCD and as(M|). 

6. Strange Spectral Function and Strange 
Quark Mass 

The spectral function for strange final states 
has been determined by ALEPH [5]: it is dom- 
inated by the vector if* (890) and higher mass 
(mostly axial-vector) resonances. The total 
rate for strange final states, using the complete 
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Figure 7. The running of as{so) obtained 
from the fit of the theoretical prediction to 
i?T-,v+A(so)below the r mass in the top plot: the 
shaded band shows the data including experi- 
mental errors, while the curves give the four-loop 
RGE evolution for two and three flavours. The 
plot below shows the evolution of the strong cou- 
pling (measured at m^) to M| predicted by QCD 
compared to the direct measurement. Flavour 
matching is accomplished at 3 loops at 2 rUc and 
2 mi, thresholds. 
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Figure 8. Differential rate for A°'^, differ- 
ence between properly normalized nonstrange 
and strange spectral functions (see text for de- 
tails). The contribution from massless pertur- 
bative QCD vanishes. To guide the eye, the 
solid line interpolates between bins of constant 
0.1 GeV^ width. 



ALEPH analyses supplemented by results from 
other experiments [33] is determined to be -B(r — > 
z/rhadrons5=_i) = (29.3 ± 1.0) 10^^^ leading to 

= 0.163 ±0.006. (11) 

Spectral moments are again useful tools to un- 
ravel the different components of the inclusive 
rate. Since we are mostly interested in the spe- 
cific contributions from the us strange final state, 
it is useful to form the difference 

\kl — ^ ukl ^ pfci ^^'->\ 

= ^^^■^=° ~ r^^^-^="^ ^'^^ 

where the flavour-independent perturbative part 
and gluon condensate cancel. Fig. 8 shows the in- 
teresting behaviour of Al^" expressed differentially 
as a function of s. The leading QCD contribution 
to A^' is a term proportional to the square of the 
strange quark mass at the t energy scale. We 
quote here the recent result from the analysis of 
Ref. [34], yielding 

m,{ml) = (120±lUp±8y„, ±19t,0 MeV(13) 
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where the dominant uncertainty is from the- 
ory, mostly because of the poor convergence be- 
haviour. 

7. Conclusions 

The decays t ^ Vr + hadrons constitute a 
clean and powerful way to study hadronic physics 
up to ^/s ~ 1.8 GeV. Beautiful resonance analy- 
ses have already been done, providing new insight 
into hadron dynamics. Probably the major sur- 
prise has been the fact that inclusive hadron pro- 
duction is well described by perturbative QCD 
with very small nonperturbative components at 
the T mass. In spite of the fact that this low- 
energy region is dominated by resonance physics, 
methods based on global quark-hadron duality 
work indeed very well. 

The ALEPH preliminary results using the full 
LEPl sample have been presented. Satisfactory 
agreement with CLEO is observed in the tttt*' de- 
cay mode. The r spectral functions have now 
reached a precision level where detailed investi- 
gations arc possible, particularly in the most in- 
teresting 7r7r° channel. The breaking of SU(2) 
symmetry can be directly determined through the 
comparison between e+e" and r spectral func- 
tions. In particular, the possibility of a p mass 
splitting between neutral and charged states is 
carefully investigated. The data tend to favour 
non-degenerate states with nip- — rUpO = (2.3 ± 
0.8) MeV. Correcting the r spectral function for 
this additional isospin breaking leaves an over- 
all normalization difference of 3.3% with e+e" 
data, thus enhancing the existing discrepancy for 
^had.LO rjnj^^ situation must be revisited with new 
high precision e+e" data [20,35]. 

The measurement of the vector and axial- 
vector spectral functions has provided the way 
for quantitative QCD analyses. These spectral 
functions are very well described in a global way 
by 0(0:^) perturbative QCD with small nonper- 
turbative components. Precise determinations of 
as agree for both spectral functions and they 
also agree with all the other determinations from 
the Z width, the rate of Z to jets and deep in- 
elastic lepton scattering. Indeed from r decays, 
as(M|)^ = 0.1202 ± 0.0027, in excellent agree- 



ment with the average from all other determina- 
tions [36], as{M^)non-T = 0.1187 ±0.0020. 

The strange spectral function yields a value 
for the strange quark mass which can be evolved 
to the usual comparison scale, ms{2 GeV) = 
(liet^O) MeV. 

The r spectral functions have been shown to 
be a privileged field for the study of QCD. 
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